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ABSTRACT 


To  further  support  the  practicability  of  the  digital  hydraulic  actuator 
as  a  control  component,  IBM  has  successfully  demonstrated  operation  of  the 
Phase  m  laboratory  model  in  a  simulated  digital  flight  control  system.  Per¬ 
formance  using  the  laboratory  model  is  compared  with  the  performance  when 
ushig  a  mathematical  model.  Recorded  transients  and  the  values  obtained  for 
a  maximum  stable  step  and  gain  change  are  presented.  The  effect  of  varying 
the  actuator  parameters  and  the  validity  of  a  linearized  model  is  also  estab¬ 
lished.  Agreement  of  the  results  shows  the  mathematical  model  to  accurately 
represent  the  actual  hardware.  Hence,  this  study  verified  the  results  of  both 
Phase  I  and  Phase  n. 
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Section  I 


INTRODUCTION 


A.  SUMMARY 


Completion  of  Phase  m  yields  results  which  further  support  the  prac¬ 
ticability  of  a  digital  hydraulic  actuator  (DHA)  for  flight  control.  A  labora¬ 
tory  model  of  the  DHA  was  operated  in  a  closed  loop  simulating  the  digital 
flight  control  of  the  Polaris  missile.  Performance  of  the  simulated  control 
system  using  the  actual  hardware  is  compared  with  the  performance  obtained 
with  the  mathematical  model.  Results  presented  show  the  mathematical 
model  to  accurately  describe  hardware  dynamics. 

Due  to  the  ram  area  ratio,  the  system  performance  could  not  be  deter¬ 
mined  with  the  required  friction  load.  However,  having  established  the  val¬ 
idity  of  the  mathematical  model,  the  performance  of  the  system  with  the 
required  ram  area  and  friction  load  was  determined  by  simulation.  The 
effect  of  varying  the  actuator  parameters  and  the  performance  for  a  "worst" 
case  parameter  condition  was  also  established.  Analog  flight  control,  which 
was  simulated  in  Phase  I,  is  used  to  provide  a  basis  for  evaluating  the  digital 
flight  control. 

A  linearized  model  of  the  DHA  is  derived.  The  model  is  shown  to  be  a 
valid  approximation  when  operation  is  not  velocity  saturated  and  quantization 
is  ignored.  Frequency  response,  phase  shift,  and  response  to  a  step  com¬ 
mand  is  used  as  a  basis  to  show  the  linearized  model  is  a  good  approximation. 

B.  PURPOSE  AND  OBJECTIVES  OF  PROGRAM 


Control  of  an  aerospace  vehicle  can  potentially  be  simplified  by  using 
the  capabilities  of  the  digital  guidance  computer  for  performing  those  control 
functions  now  being  accomplished  by  analog  techniques.  To  study  the  impact 
of  digitally  controlled  actuators  on  flight  control  systems  and  to  ascertain 
those  requirements  of  the  DHA,  ONR  initiated  Contract  No.  Nonr  4375(00). 
The  total  objective  of  this  program  is  to  demonstrate  the  feasibility  of  a  DHA 
as  a  flight  control  component.  The  program  consists  of  three  phases  whose 
objectives  are: 

•  Phase  I  -  To  ascertain  those  DHA  requirements  that  are  necessary 
in  an  aerospace  application  such  as  Polaris 


•  Phase  n  -  To  design  and  evaluate  DHA  meeting  the  requirements 
found  necessary  in  Phase  I 

e  Phase  in  -  To  demonstrate  the  performance  of  the  DHA  hardware 
in  a  simulated  digital  flight  control  system  and  to  establish  the 
accuracy  of  the  mathematical  model. 

The  use  of  an  incremental  DHA  will  result  in: 

•  System  simplification  by  the  elimination  of  digital- to- analog  con¬ 
version  equipment  required  for  communicating  between  the  digital 
computer  and  the  actuator 

e  Elimination  of  an  output  transducer  and  a  summing  amplifier  in  the 
analog  hydraulic  servo  loop 

e  Substitution  of  bi-  stable  control  valves  for  the  analog  control  valves, 
the  latter  being  more  susceptible  to  null  shift  and  oil  contamination 

e  Increased  versatility  to  perform  such  operations  as  programmed 
gain  changes  or  adaptive  functions 

•  A  significant  step  toward  realizing  an  all-digital  control  system. 

C.  DESCRIPTION 


The  IBM  Digital  Metering  Unit  complemented  by  a  hydraulic  actuator 
make  up  the  DHA.  The  metering  unit  actuated  by  digital  commands  produces 
a  required  change  of  nosition  by  the  transfer  of  calibrated  fluid  volumes. 

Each  digital  command  produces  a  volume  of  fluid  equal  to  one  quantum  of  po¬ 
sition.  The  DHA  thus  provides  a  capability  of  transforming  digital  commands 
without  conversion  into  related  output  positions.  A  detailed  description  of  the 
operational  concept  is  presented  in  Appendix  A. 

A  mathematical  model  of  the  DHA  is  compared  with  the  actual  hardware 
by  the  arrangement  depicted  schematically  in  Figure  1-1.  The  arrangement 
provides  a  capability  of  monitoring  the  unused  model.  Equations  describing 
the  flight  control  system  and  the  DHA  are  implemented  by  a  hybrid  combina¬ 
tion  of  an  IBM  709  EDPM  and  a  EAI  Analog  Computer.  The  arrangement  of 
the  hardware  is  shown  in  Figure  1-2. 


Figure  1-2.  Laboratory  Arrangement 


D.  CONCLUSIONS 


Completion  of  Phase  m  yielded  the  following  conclusions: 

•  The  mathematical  model  of  the  DHA  is  an  accurate  representation 
of  the  actual  hardware 

•  Results  obtained  by  the  Phase  I  simulation  are  valid 

•  Parameter  values  to  attain  a  required  system  performance  can  be 
accurately  determined  with  the  mathematical  model  is  a  simulated 
system. 

•  A  DHA  is  feasible  for  an  aerospace  application,  such  as  Polaris 

•  For  an  application  such  as  the  Polaris,  the  acceptable  tolerance 
on  the  leakage  parameters,  qL,  and  q^  ,  is  relatively  large. 
Hence  leakage  variation  due  to  production  tolerances  and  wear 
should  not  present  a  problem 

•  The  suitability  of  the  DHA  to  a  given  application  can  be  initially 
determined  by  means  of  a  linearized  model. 

E.  RECOMMENDATIONS 


The  DHA  for  this  study  was  tailored  to  the  requirements  of  the  Polaris 
Missile.  Other  applications  may  require  a  DHA  capable  of  producing  a 
smaller  resolution  and  a  higher  velocity.  For  some  applications  such  as 
aircraft,  a  two- speed  actuator  may  be  necessary  to  achieve  both  high  reso¬ 
lution  and  high  velocity  of  the  controlled  surface. 

The  following  is  recommended  for  future  study: 

•  Ascertain  by  simulation  the  requirements  and  performance  when 
using  a  DHA  for  the  control  of  an  aircraft,  including  the  effect  of 
the  oscillatory  motion  of  the  control  surface  during  steady  state 
on  the  vibratory  modes  of  the  structure. 

•  Evaluate  the  performance  of  a  two- speed  actuator  having  the  capa¬ 
bility  of  changing  the  direction  of  drift* 


♦Raider,  Scott,  "Digital  Hydraulic  Study  Engineering  Summary  Report  Phase 
H;  Contract  No.  Nonr  4375(00),  IBM  No.  66-928-7. 


•  Evaluate  the  design  of  the  laboratory  model  of  the  DHA.  Based  on 
the  results  design  a  DHA  which  is  representative  of  a  deliverable 
item,  meeting  reliability  and  environment  specifications. 

•  Design  of  a  metering  unit  which  uses  a  double-ended  ram.  This 
will  increase  the  output  power  of  the  actuator  for  any  selected 
flow  and  pressure. 
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Section  n 


TECHNICAL  APPROACH 


Phis  section  describes  the  equipment,  and  implementation  of  the  equa¬ 
tions.  The  Phase  m  Study  consisted  of  the  following  effort: 

•  Design  logic  interface  between  the  hardware  and  simulated 
flight  control  system 

•  Implement  equations  defining  the  DHA  and  flight  control  system 

•  Compare  performance  of  simulated  control  system  when  using 
actual  hardware  with  the  performance  when  using  the  mathe¬ 
matical  model 

•  Establish  the  performance  of  the  Polaris  missile  when  using  a 
DHA  which  provides  a  force  independent  of  the  direction  of  mo¬ 
tion  and  opposed  by  a  600-pound  coulomb  friction  load 

•  Ascertain  the  effect  of  parameters  on  performance  of  the  simu¬ 
lated  system 

•  Determine  the  validity  of  approximating  a  DHA  and  digital  con¬ 
trol  system  by  a  linearized  model. 

Phase  HI  is  consistent  with  Phase  I  in  that: 

•  Missile  dynamics  are  when  the  missile  is  subjected  to  maximum 
aerodynamic  pressure  which  corresponds  to  the  least  controllable 
portion  of  the  flight  trajectory 

•  The  yaw  angle,  the  loop  of  highest  gain,  is  studied 

•  Feedback  components  are  assumed  to  be  ideal  in  that  no  error 
is  contributed  to  the  system  by  the  resolution  or  drift  of  the 
attitude  sensors 


•  The  transfer  function  of  the  Polaris  missile  is  derived  for  the  short 
period  oscillation  mode,  assuming: 

—  Rigid  Structure 

—  Stationary  environment 

—  Small  perturbations 

—  Missile  aligned  with  the  equilibrium  position  of  the  velocity 
vector  (zero  sideslip) 

—  Constant  mass. 

A.  EQUIPMENT  -  ARRANGEMENT  AND  DESC  RIPTION 

Arrangement  of  the  equipment  for  the  Phase  m  simulation  is  shown  in 
Figure  2- 1.  Two  modes  of  operation  are  provided,  the  analog  mode  and  sim¬ 
ulation  mode.  The  analog  mode  permits  the  DHA  to  be  operationally  checked 
without  the  computer.  A  commanded  position  is  compared  with  the  actual  po¬ 
sition  to  produce  an  error  signal  which  is  converted  into  digital  signals;  a 
sign,  and  a  volume  command.  The  simulation  mode  permits  the  DHA  to  be 
operated  open  loop  to  control  a  simulated  plant  (the  Polaris  missile).  Po¬ 
sition  error  is  determined  by  the  control  loop  about  the  plant. 

1.  VALVE  CONTROL  LOGIC 

The  control  logic  permits  evaluation,  in  either  the  simulation  or  analog 
mode,  of  a  single- speed  actuator,  two- speed  actuator  or  an  actuator  for 
achieving  accurate  positioning,  i.  e. ,  operation  at  the  discontinuous  point  of 
the  flow  curve  (Figure  2-  2).  Schematic  and  detailed  logic  is  presented  in 
Appendix  B. 

Logic  required  for  the  Phase  HI  Study  is  depicted  ir.  Figure  2-3.  The 
logic  requires  that  a  sign  command  be  received  in  parallel  with  each  volume 
command.  Sign  commands  following  the  last  sign  change  produce  Uie  gating 
signal,  A.  A  volume  command  when  accompanied  with  the  signal,  A,  is  gated 
to  the  actuating  driver.  The  actuating  driver,  upon  receipt  of  a  digital  com¬ 
mand,  de-energizes  the  valve  latch  to  initiate  the  transfer  of  a  calibrated 
fluid  volume.  When  changing  the  direction  of  motion,  the  signal, 5,  does  not 
occur,  and  the  volume  command  in  parallel  with  the  sign  command  is  inhib¬ 
ited.  The  sign  command  switches  the  sign  driver,  and  the  transfer  of  a  cal¬ 
ibrated  volume  follows  the  switching  of  the  sign  valve.  Sign  commands  follow¬ 
ing  the  last  sign  change  are  transmitted  to  the  same  sign  driver.  However, 
the  sign  valve  is  now  controlled  by  the  sign  driver  associated  with  the  latch 
magnets  holding  the  armature. 


Simulation 


Figure  2-1.  Arrangement  of  Equipment 


InUt  Fluid  Temperature  —  100#F  9 


Figure  2-2.  Flow  Curve  for  Digital  Unit 

2.  POSITION  CONTROL  LOGIC 

The  position  control  logic  shown  in  Figure  2-3  is  used  when  operating 
in  the  analog  mode  to  produce  the  sign  and  volume  commands.  The  error 
detector  determines  the  correction  required  by  predicting  the  error  that  will 
exist  after  effecting  the  last  command.  This  is  accomplished  by  updating 
the  error  signal,  based  on  present  position,  with  the  signal  produced  by  the 
current  driver,  representing  the  position  change  commanded.  Hence,  the 
time-lag  due  to  the  time  required  to  position  the  load  does  not  introduce  a 
false  error.  The  error  detector  outputs  (sign  and  volume  commands)  are 
gated  to  the  pulse  generators  to  produce  pulses  compatible  with  the  valve 
logic  input  requirements. 

3.  DIGITAL  HYDRAULIC  ACTUATOR 

As  previously  stated,  the  DHA  consists  of  the  prototype  digital  metering 
unit  and  a  commercial  actuator  cylinder.  The  following  modifications  were 
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Figure  2-  3.  Phase  m  Control  Logic 


incorporated  to  improve  performance  and  overcome  operation  problems  ex¬ 
perienced  during  Phase  II: 

•  Relative  motion  between  the  pilot  spools  and  armature  bearing 
play  was  minimized  by  the  use  of  shims  and  removing  the  bearing 
located  on  the  eccentric 

•  New  magnets  were  designed  to  provide  a  larger  holding  force  and 
to  eliminate  the  problem  of  lead  wire  breakage 

•  A  new  manifold  was  designed  to  eliminate  the  excessive  leakage 
due  to  the  location  of  the  port  to  ram.  In  addition,  the  manifold 
provides  ports  for  an  additional  valve  unit  and  calibrating  cylinder. 
Hence  a  capability  for  two- speed  operation. 

A  commercial  cylinder  having  a  ram  area  of  0.  994  sauare  inch  and 
an  area  ratio  of  1.  59:1  was  used.  Since  this  ratio  differs  from  that  re¬ 
quired,  the  force  produced  is  not  the  same  for  positive  and  negative  motion. 

4.  COULOMB  FRICTION 

Coulomb  friction  is  generated  by  the  bearing  and  seals  when  firing  the 
nozzle.  To  simulate  coulomb  friction,  a  phenolic  pad  was  hydraulically 
loaded  against  an  aluminum  drum  segment  connected  to  the  nozzle.  Radial 
loading  of  the  bearings  was  minimized  by  placing  the  back  of  the  drum  segment 
against  the  outer  race. 

Because  the  area  ratio  of  the  ram  limited  the  force  capability,  the 
required  coulomb  friction  could  not  be  used  when  operating  with  the  hardware. 

B.  METHOD  OF  SIMULATION 

The  mathematical  model  for  the  Phase  in  Study  of  the  digital  autopilot, 
DHA  and  missile  dynamics  is  shown  in  Figure  2-  4. 

1.  DIGITAL  AUTOPILOT 

The  digital  autopilot  is  programmed  as  shown  in  Figure  2-5.  The  yaw 
command,  Ny(t),  can  be  either  a  step  |Ny  |  at  t  =  0  or  j  Ny|  sin  wt. 


Figure  2-4.  Mathematical  Model  of  Digital  Autopilot. 
Digital  Actuator,  and  Missile  Dynamics 
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The  pulse  logic  box,  which  is  depicted  in  detail  in  Figure  2-  6,  has  been 
generalized  to  allow  for  the  following  autopilot  options: 

•  Single- Speed  Operation  —  Single- speed  operation  was  used  for  the 
Phase  in  Study.  This  option  is  attained  by  setting  eg  =  e^.  Note 
that  the  condition  imposed  on  c  to  select  a  small  volume  command 
is  then  e^  >|c  |  >  e8  which  cannot  be  satisfied  since  es  =  e^. 

•  Two- Speed  Operation  —  This  option  makes  it  possible  to  simulate 
an  actuator  having  two  volume  sizes,  a  large  volume  for  slew,  and 
a  small  volume  for  fine  quantization. 
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Figure  2-6.  Pulse  Logic 


•  Accurate  Positioning  (Precision  Mode)  —  This  option  provides  the 
capability  of  changing  the  direction  of  ram  motion  without  displacing 
a  calibrated  volume. 

If  |c  |  <  efl  the  drift  option  allows  the  drift  polarity  to  be  de¬ 
termined  by  the  polarity  of  "Input”,  <  ,  or  the  last  command. 

•  Feedback  Flexibility  —  The  polarity  of  c  may  be  the  opposite  of  the 
DHA  sign  pulse  for  |  cl  <  e8.  To  provide  the  flexibility  of  attaining 
a  feedback  (FB)  polar  lily  equal  to  the  polarity  of  cor  the  polarity  of 
the  sign  pulse,  FB  is  set  equal  to 

8n  sgn  c  +  8p  sgn  (sign  pulse). 

For  small  volume  commands 

FB  =  8  s  sgn  c  , 

and  for  large  volume  commands 
FB  =  8b  sgn  c. 

Hammering  of  the  actuator  into  its  stops  by  the  control  computer  is 
prevented  during  check  out  and  stability  check  runs  by  inhibiting  the  volume 
commands  if  the  ram  position  magnitude,  |X]J  ,  exceeds  it-  A  value 
of  XLim  was  selected  so  that  |Xi  |  =  Xnm  when  the  ram  is  0.  3  inch  from  the 
stop.  (This  distance  allows  for  the  displacement  produced  by  the  volume  stored  by 
the  auxiliary  pistons. )  Since  the  plant  may  recover  from  short  periods  of  8a 
saturation,  two  limit  options  are  provided.  Limit  option  1  assumes  the  con¬ 
trol  computer  has  no  feedback  from  the  actuator  and  does  not  alter  FB.  Op¬ 
tion  2  allows  for  a  special  FB  equal  to  8l  sgnc  assuming  some  sort  of  limit 
detection.  In  summary,  Phase  I  is  duplicated  if: 

e.  =  e  =  1  quantum  +  hysteresis  _  q  ggj  for  a  quanti2er  hysteresis 

b  8  2  of  8_  /yl 


drift  option  =  3 
limit  option  =  1 


'  m 


2.  DIGITAL  HYDRAULIC  ACTUATOR  (DIGITAL  PORTION) 

After  transmitting  the  DHA  pulses,  the  digital  computer  calculates 
Qpl,  Qp2>  Qp3  and  AV  as  shown  in  Figure  2-7.  Where: 

Qpl  is  the  flow  for  t  in  the  interval  (tv ,  tfS ) 

Qp2  is  the  flow  for  t  in  the  interval  (t{8 ,  tfg ) 

Qp3  is  the  flow  for  t  in  the  interval  (tf  3 ,  ty ) 

AV  =  ty  +  T 

/  V' 

tv  (i  ■  1,  2,  3) 


Figure  2-7.  Flow  Calculations 


With  reference  to  the  timing  chart  shown  in  Figure  2-8,  the  Qp. 

(i  =  1,  2,  3)  are  converted  from  digital  to  analog  at  the  proper  time  and  fed 
to  an  analog  integrator  to  get  Vp.  To  prevent  drift  in  this  open- loop  integra¬ 
tion,  the  AVs  are  summed  to  form  a  digital  representation  of  Vp.  A  drift 
correction  (Qcorr)  is  added  to  the  Qpi  to  force  the  analog  vahrn  of  Vp  to  match 
the  digital  value  Vdig,  at  sampling  times  (Figure  2-9). 

Transfer  of  the  volume  commanded  during  the  interval  T(n  _  is  com¬ 
pleted  during  the  time  interval  T/n\ .  Therefore,  to  make  VdJg  assume  values 
correct  at  the  analog- to- digital  version  time,  T^/p,  of  Vp,AV  is  divided  so 
that 


AV 


t  +  T,  x 
.v  (n) 

/  Qpdt  =  AV1+AV2 


where 


/ 


A/D 

Q 


+  T 
dt 


AV2 


ty  +  T 

/«Pdt 

WT 


Thus  during  the  nth  DHA  cycle,  the  volume  added  to  the  integration 
summer  is: 


AV'  =  AVj  (n)  +  AV2  (n  -  1) 


or 


AV’  = 


Vd+t 

/  <V‘ 

lA/D 
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DHA  Period  T 


Figure  2-8.  Timing  Chart 


Figure  2-  9.  Digital  Integration 


ure  2 


The  limiting  values,  V^ay  and  Vm^,  for  the  digital  integrator  (Fig- 
:-9)  are  variable  and  are  obtained  from  the  following  relationships: 


which  is  depicted  by  block  diagram  in  Figure  2-  10.  Since  Vj  has  a  fixed 
limit,  VLim> 

-VLim  <  Vt  K  VLim 
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Figure  2-10.  Method  of  Obtaining  Limiting  Values  for  the 
Digital  Integrator 
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Therefore,  Vs,  an  analog  quantity,  is  converted  from  analog  to  digital  form, 
and  the  limits  on  V^g  are  computed  by  the  digital  computer  as 

v  =  -V  +  V 
vmin  vLim  VB 


max 


VLim  +  V 


3.  DIGITAL  HYDRAULIC  ACTUATOR  (ANALOG  PORTION) 

The  analog  simulation  of  the  actuator  dynamics  is  programmed  as  in 
Phase  I,  (Figure  2-11)  with  the  exception  of  the  coulomb  friction.  As  a  re¬ 
sult  of  simulating  in  real  time  (Phase  I  simulation  was  1/10  as  fast)  the 
limited  frequency  response,  and  amplifier  phase  shift  caused  the  coulomb 
feedback  loop  to  oscillate  when  used  as  a  simple  bang- bang  amplifier  as  ill¬ 
ustrated  in  Figure  2-12.  Hence  during  a  steady  state,  large- amplitude  os¬ 
cillations  of  the  pressure,  resulted.  To  overcome  this  problem,  a 
latching  electronic  swtich  is  used  at  the  input  of  the  X  integrator  (Figure  2- 
13). 


Figure  2-13.  Simulation  of  Coulomb  Friction  (Phase  HI) 


The  logic  that  operates  switch  S  is: 

•  S  closes  and  latches  whenever  |  P^AL  |  >F 

•  S  unlatches  when  both  j  PdAL  j  <  F 


X  =  0 


This  prevents  AF  from  acting  as  an  accelerating  force  because  of 
phase  shift. 

4.  COMPARISON 

This  method  of  simulating  the  digital  autopilot  and  DHA  for  Phase  IQ 
differed  from  Phase  I.  To  establish  that  the  results  obtained  with  the  present 
model  agreed  with  the  results  obtained  during  Phase  I,  several  simulation 
runs  depicting  significant  characteristics  were  rerun.  Parameter  values  used 
for  the  reruns  were  identical  to  those  used  for  the  Phase  I  Study  except  for 
the  fluid  flow  time,  Tfg .  The  parameter  Tfg  is  not  readily  modified  in  the 
program;  therefore,  Tfg  was  programmed  fo**  the  value  obtained  by  test  and 
does  not  significantly  influence  the  results. 

The  results  attained  are  compared  in  Table  2- 1  with  the  results  of 
Phase  I.  The  agreement  of  the  two  methods  A  implementing  the  equations 
is  further  illustrated  by  the  recorded  transients  shown  in  Figure  2- 14. 
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Time  (seconds) 
a)  Phase  I  Simulation 


Time  (seconds) 
b)  Phase  III  Simulation 


Figure  2-14.  Recorded  Transients  for  \[/  =  -1/2 
(Run  207  Table  2-1) 
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Section  m 


MATHEMATICAL  MODEL  EVALUATION 


Performance  of  the  simulated  control  loop  when  using  the  mathematical 
model  of  the  DHA  is  compared  with  the  performance  when  using  the  labora¬ 
tory  hardware.  Accuracy  of  the  model  is  established  by  correspondence  of 
the  recorded  ram  motion  and  agreement  of  the  results  for: 


•  Maximum  stable  step 

•  Response  to  a  step  command,  N^  =  10 

•  Response  with  an  initial  yaw  condition 

•  Gain  margin  for  initial  conditions  less  than  one  quantum 

•  Limit  cycle 

•  Frequency  response  and  phase  shift. 


Plant  transients  presented  in  this  section  are  of  significance  only  for 
the  purpose  of  comparison  since  the  ram  used  did  not  produce  the  same  force 
for  both  directions  of  motion. 

With  the  DHA  in  the  loop  the  initial  condition  of  the  ram  is  obtained  by 
allowing  the  ram  to  drift  from  a  positive  to  a  negative  value.  The  position 
is  sampled  at  the  operating  frequency  (82. 5  Hz)  and  the  run  is  started  when 
a  zero  or  negative  ram  position  is  detected. 


The  unpredictable  behavior  of  an  unstable  plant  with  zero  initial  con¬ 
ditions,  quantization  effects,  and  the  method  of  initializing  the  position  of 
the  physical  load  will  produce  variations  in  the  plant  response.  Therefore, 
small  /ariations  of  the  plant  transient  when  using  the  mathematical  model 
from  the  transient  generated  when  using  the  hardware  does  not  necessarily 
indicate  an  inaccuracy  of  the  mathematical  model. 

A.  PARAMETER  VALUES 

Parameter  values  defining  the  actuator  and  load  dynamics  are  shown 
in  Table  3-1.  The  parameter  values  were  established  by  test  when  practi¬ 
cable.  Values  not  measured  are  indicated.  Accurate  values  are  not  readily 


Table  3-1 


PARAMETER  VALUES 


Parameter 

Units 

Value 

Parameter 

Units 

Value 

AF 

lb 

380 

®b 

deg 

0.681 

M 

lb-s2/in 

0.25 

Quant.  Hys. 

deg 

Sq/4 

R 

in 

8.125 

8b 

deg 

1.09 

*0 

T 

lb-in“5 

26,600 

K 

— 

2"7 

*c 

lb-s/'in 

12 

in3/s 

0. 107 

qLC 

in3/s 

-0.006 

^LC* 

in3/s 

-0. 0165 

VCB 

In3 

0.152 

*C 

Hz 

82.5 

TV 

in'3 

5.5 

lf 

In"3 

10 

PM 

lbf/in2 

+540 

-960 

PR 

lbf/in2 

130 

Ph 

lbf/in2 

2,130 

.  2 
in 

0.  994 

al/ar 

— 

1.59 

♦Value  Approximated 


B.  RESULTS 


Simulation  results  when  using  the  parameter  values  of  Table  3-1  are 
shown  in  Table  3-2.  Table  3-2  also  shows  the  results  when  changing  the 
parameter  value  of  Table  3-1  to  that  indicated. 

Table  3-2 

PERFORMANCE  OF  SIMULATED  FLIGHT  CONTROL  SYSTEM 


Parameter 

Value 

Maximum 

Stable 

— 

Unstable 

Step 

ii 

£ 

o' 

ii 

O  ° 
US 

Reference: 
Figure  No. 

Step 

Stable 

Unstable 

Ny  =  10 

Limit 

Cycle 

Hardware 

Digital 

Hydraulic 

Actuator 

11 

-8 

12 

-9 

4 

4.5 

3-2 

3-8 

3-5 

Mathematical 

Model 

Table  3-1 

11 

-8 

12 

-9 

4 

4.5 

3-6 

3-7 

Mathematical 
Model  with: 

C  =  6 

11 

-9 

12 

-10 

4 

4.5 

3-9 

3-11  a 

C  =  18 

11 

-6 

12 

-7 

3.5 

4 

3-10 

3-11  b 

Y  =  13, 300 

11 

-8 

12 

-9 

4 

4.5 

3-12  a 

3-13  a 

^  =  100, 000 

11 

-8 

12 

-9 

4 

4.5 

3-12  b 

3-13  b 

o 

II 

< 

11 

-9 

12 

-10 

4.5 

5 

Recorded  transients  which  further  support  the  accuracy  of  the  mathema¬ 
tical  model  and  the  effect  when  varying  the  parameters  are  shown. 

As  shown  by  Table  3-1,  the  maximum  stable  step  is  dependent  on  the 
sign  of  the  steering  pulse  command  Ny.  The  difference  between  the  maxi¬ 
mum  positive  and  maximum  negative  stable  step  is  a  consequence  of  the  ram 
force  being  different  for  positive  and  negative  motion.  The  higher  negative 
force  produces  a  larger  effective  moment  due  to  faster  response  and  in¬ 
creased  overshoot.  Recorded  transients  of  the  ram  for  positive  and  negative 
yaw  step  commands  are  shown  in  Figure  3-1.  For  the  negative  step  command 
the  ram  is  seen  to  be  velocity  saturated  during  the  positive  motion. 

Figure  3-2  shows  the  recorded  ram  motion  of  the  laboratory  model 
and  mathematical  model  for  a  yaw  step  command  when  closing  the  simulated 
loop  with  the  laboratory  model  (Figure  3-3).  It  is  significant  that  the  tran¬ 
sients  are  parallel.  With  the  mathematical  model  operated  open  loop,  and 
error  in  the  leakage  rate,  qLC»  cause  the  ram  transient  for  the  mathe¬ 
matical  model  to  converge  or  diverge  relative  to  the  transient  recorded  for 
the  actual  ram.  Hence  the  parallel  transients  verify  the  value  of  to  be 
reasonably  correct. 

Although  the  agreement  is  good  the  transients  differ  when  changing 
from  negative  ram  motion  to  positive  ram  motion,  detail  A,  Figure  3-2. 

This  difference  is  due  to  the  velocity  produced  with  the  mathematical  model 
for  negative  motion  being  greater  than  that  achieved  with  the  hardware.  A 
lumped  parameter  approximation  is  used  to  describe  damping  in  the  mathe¬ 
matical  model.  Accuracy  of  the  approximation  is  decreased  when  the  damp¬ 
ing  is  significantly  different  for  positive  and  negative  motion.  With  the 
larger  accelerating  force  for  negative  commands,  Figure  3-4,  a  larger  force 
acts  on  the  ram  before  the  auxiliary  piston  is  displaced.  Until  the  calibrating 
piston  stroke  is  completed,  the  return  lines  effect  damping  of  the  ram  motion. 
Hence,  because  of  the  unequal  ram  force,  the  damping  realized  when  initiating 
negative  motion  is  greate  than  that  produced  by  the  mathematical  model. 

The  value  used  for  the  damping  coefficient  (C  =  12)  better  represents 
the  damping  for  positive  motion.  This  is  shown  by  the  agreement  when  the 
order  of  ram  motion  is  reversed,  see  detail  B,  Figure  3-2. 

Limit  cycles  produced  with  the  hardware  and  mathematical  model  in 
parallel  are  shown  in  Figure  3-5.  With  the  simulated  loop  closed  about  the 
hardware,  the  higher  velocity  achieved  with  the  mathematical  model  for 
negative  commands  can  be  seen.  Also  of  interest  is  the  "stepping”  charac¬ 
teristic  produced  with  the  mathematical  model  during  drift  (see  detail  A  of 
Figure  3-5).  This  characteristic  is  due  to  the  method  of  implementing 
coulomb  friction.  Steps  occur  when  the  condition  AL  £  AF  is  obtained, 
thus  closing  the  switch  (Figure  2-13).  Oscillations  of  Pd  which  result 
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Figure  3-1.  Recorded  Motion  of  Ram  with  Digital  Hydraulic 
Actuator  in  Simulated  Control  Loop  (Af  =  60  lbs) 


Time  (seconds) 


Figure  3-2.  Ram  Motion  for  a  Yaw  Step  Command  (Ny  =  10) 

create  the  conditions  required  to  open  the  switch,  i.  e. ,  Pd  Al<  AF,  Xi  =  0. 
The  stepping  characteristic  can  also  be  noted  in  Figure  3-13  showing  the 
effect  of  the  parameter  (3/V  on  the  limit  cycle.  With  the  increased  resonance 
frequency  (larger  /3/V)  the  stepping  frequency  increases. 

The  ram  transient  for  a  yaw  step  command  using  the  mathematical 
model  in  the  simulated  loop  is  shown  in  Figure  3-6.  Limit  cycle  characteris¬ 
tics  are  shown  in  Figure  3-7.  Comparison  with  Figures  3-2  and  3-5  shows 
the  recorded  ram  motion  to  exhibit  the  same  characteristics  as  that  produced 
by  the  hardware.  Comparison  of  the  plant  transients  shown  in  Figures  3-6 
with  Figure  3-8  show  the  response  characteristics  of  the  simulated  loop  with 
the  mathematical  model  of  the  DHA  to  be  the  same  as  that  achieved  with  the 
laboratory  hardware. 
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Figure  3-3.  Simulation  Arrangement 


Figure  3-4.  Block  Diagram  for  DHA  Load  Dynamics 
1.  EFFECT  OF  THE  PARAMETER  C 

The  effect  of  parameter  C  on  loop  performance  is  shown  by  the  results 
in  Table  3-2.  The  values  of  C  represent  a  variation  ±50  percent  from  that 
used  when  comparing  the  mathematical  model  with  the  hardware.  For  both 
increased  and  decreased  damping,  the  maximum  negative  stable  step  differed 
from  the  result  obtained  with  the  laboratory  model.  However,  the  maximum 
positive  stable  step  was  not  changed.  As  shown  by  Figure  3-1  negative  ram 
motion  is  required  to  stabilize  the  missile  after  initiating  a  positive  step 
command.  For  the  values  of  damping  investigated  the  relation: 

Vc  .  PmAL-AF 

- C - 

existed  with  negative  motion.  Hence  the  slew  velocity  required  to  attain  the 
maximum  positive  stable  step  was  not  effected  by  the  values  of  C  investigatec 


Figure  3-5.  Recorded  Transients  Showing  Actual  Motion  of  Digital 
Hydraulic  Actuator  and  the  Motion  Generated  with  the  Mathematical  Model 


Figures  3-9  and  3-10  show  the  response  to  a  yaw  step  command  for  C  =  6  and 
C  =  18.  The  variation  of  the  damping  parameter  produces  no  significant 
change  in  the  transient  response.  The  effect  of  the  damping  on  the  ram 
,  motion  is  shown  in  Figure  3-11.  As  to  be  expected,  the  overshoot  of  ram 

|  position  is  greater  for  decreased  damping.  However,  the  limit  cycle 

I  characteristics  as  shown  in  Figures  3-9  and  3-10  do  not  significantly  differ. 

j  2.  EFFECT  OF  THE  PARAMETER  /3/V 

The  effect  of  the  parameter  £/V  was  determined  for  values  of  13, 300 
lb/ in®  and  100, 000  lb/ in®.  As  indicated  by  Table  3-2  the  variations  had  no 
effect  on  gain  margin  or  maximum  stable  step.  Transients  ret  srded  for  a 
yaw  command  Ny  =  10  are  shown  in  Figure  3-12  for  the  variations  of  the 
parameter  /3  /V.  Both  transients  are  in  satisfactory  agreement  with  the 
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Figure  3-6.  Recorded  Transients  for  a  Yaw  Step  Command  (Ny  =  10) 

Using  Mathematical  Model 
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Figure  3-7.  Limit  Cycle  Generated  with  the  Mathematical  Model  in 
the  Simulated  Control  System 
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Figure  3-8.  Recorded  Transients  for  a  Yaw  Step  Command  (Ny  -  10) 
Using  Actual  Hardware  with  Mathematical  Model  Operated  Open  Loop 

(Ref  -  Figure  3-3) 
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Figure  3-9.  Response  to  a  Yaw  Step  Command 
with  C  =  6  lb-s/ in 

transients  recorded  for  the  hardware  and  the  mathematical  model.  The 
difference  in  the  settling  time  of  the  plant  relative  to  that  depicted  in  Figure 
3-6  is  noted.  However,  due  to  the  correspondence  between  the  transients 
of  Figure  3-12  the  difference  in  settling  time  is  attributed  to  a  change  of 
conditions  within  the  simulated  loop. 

With  reference  to  Figure  3-13,  the  effect  of  the  parameter  )9/V  on 
drift  characteristics  of  the  ram  and  ram  overshoot  is  shown.  Comparison 
with  Figure  3-7  shows  the  value  for  0/V  =  26, 600  lb/in^  to  provide  better 
agreement  with  ihe  hardware. 
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Figure  3-10.  Response  to  a  Yaw  Step  Command 
with  C  =  18  lb-s/in 

The  parameter  &/V  for  the  values  stated  produced  little  change  in 
the  steady  state  frequency  of  the  ram,  see  Figure  3-13. 

3.  GAIN-PHASE 

Closed  loop  frequency  response  and  phase  characteristics  of  the 
simulated  loop  were  obtained  when  using  the  laboratory  model  and  the 
hardware.  Results  are  depicted  in  Figures  3-14  and  3-15.  Inputs  were 
restricted  in  amplitude  to  eliminate  the  nonlinear  effect  of  velocity  saturation. 
Differences  when  comparing  results  are  attributed  more  to  the  accuracy  of 
measuring  the  data  than  to  a  difference  in  the  models. 

The  agreement  of  results  further  establishes  the  accuracy  of  the 
mathematical  model  describing  the  digital  hydraulic  actuator. 
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Figure  3-11.  Limit  Cycle 
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Figure  3-15.  Closed  Loop  Amplitude  and  Phase  Characteristics 
of  Simulated  Digital  Flight  Control  System  Using 
Mathematical  Model  (See  Table  3-1) 


Section  IV 


DIGITAL  FLIGHT  CONTROL  SIMULATION 


A.  DESCRIPTION 

Performance  of  the  Polaris  missile  was  determined  when  using  a  DHA 
which  provides  the  same  force  for  both  directions  of  motion.  A  nozzle  torque 
independent  of  velocity,  producing  an  axial  load  of  600  pounds,  is  assumed. 
Results  were  obtained  with  the  mathematical  model  for  the  actuator  param¬ 
eters  of  the  actual  hardware  and  a  ram  area  ratio  of  1. 885:1.  The  equivalent 
mass  load  differs  from  that  used  when  evaluating  the  mathematical  model 
due  to  the  longer  moment  arm  used  with  the  hardware.  Parameter  values 
are  summarized  in  Table  4-1. 

The  effect  of  the  leakage  parameters  and  the  gain,  K,  of  the  autopilot 
compensation  is  shown.  Performance  is  evaluated  by  comparison  with  the 
performance  of  the  analog  flight  control  simulated  during  Phase  I. 

B.  RESULTS 


Simulation  results  for  the  parameter  values  of  Table  4-1  are  shown  in 
Table  4-2.  Also  shown  are  results  for  the  parameter  values  indicated. 

1.  EFFECT  OF  LEAKAGE  PARAMETERS 

As  shown  by  the  results  of  Table  4-2  system  performance  of  the  flight 
control  system  was  not  decreased  for  the  assumed  values  of  the  leakage 
parameters,  qL,  and  qy  q.  A  worst-case  condition  was  assumed  by  using 
the  larger  of  the  values  tor  both  q.  and  q^-,.  Values  greater  than  shown 
could  be  used,  however,  the  variation  from  that  of  the  actual  hardware  is 
believed  to  represent  a  reasonable  tolerance.  Figure  4-1  shows  the  recorded 
transients  for  a  yaw  step  command  using  the  parameter  values  of  Table  4-1. 
The  transient  of  Xj  in  Figure  4-1,  is  typical  of  that  recorded  for  the  leakage 
conditions  shown  in  Table  4-2. 

The  recorded  transients  for  an  initial  condition  0  -  - 1/2  degree  are 
shown  in  Figure  4-2  for  the  parameter  values  of  Table  4-1.  Results  for  a  no 
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Table  4-1 


PARAMETER  VALUES  FOR  SIMULATION  OF 
DIGITAL  FLIGHT  CONTROL 


Parameter 

Units 

Value 

Parameter 

Units 

Value 

Af 

lb 

600 

eb 

deg 

0. 681 

M 

lb-s2/in 

0. 296 

Quant.  Hys. 

deg 

Sq/4 

R 

in 

8.00 

8b 

deg 

1.09 

0 

V 

lb-  in"  3 

26,600 

K 

— 

2-7 

C 

lb- s/in 

12 

«L 

in3/s 

0.107 

*LC 

in3/s 

-0.006 

^LC 

in3/s 

-0.0165 

VCB 

in3 

0.152 

fc 

Hz 

82.5 

Tv 

in'3 

5.5 

Tf 

in"3 

10 

PM 

psi 

+750 

PR 

psi 

130 

Ph 

psi 

2130 

al 

in2 

1.00 

al/ar 

— 

1.885 

-2 


leakage  condition  are  shown  in  Figure  4-3  and  for  the  worst-case  leakage  in 
Figure  4-4.  Power  during  steady  state  can  be  expressed  by: 


W 


X1  al  P' 


,  ar 


c’ 


where  f  is  the  actuator  frequency. 

Table  4-2 

SUMMARY  OF  RESULTS  - 
SIMULATION  OF  DIGITAL  FLIGHT  CONTROL 


Parameter 

Maximum 

Stable 

Step 

Unstable 

Stable 
.  *G 

(Ny  = 

Unstable 

Reference  Figure  No. 

Value 

Step 

Ny  =  10 

KSE 

Table  4-1 

10 

11 

4.0 

■1 

4-1 

4-2 

qL  ■  • 150 

10 

11 

4.5 

D 

qLC  =  0165 

+11 

-10 

12 

-11 

4.5 

5.0 

’l'Ilc  =  0 

10 

11 

4.5 

■a 

4-3 

qT  =  .  150  1 

11 

12 

1 

qLC=-0165J 

-10 

-11 

4.5 

B 

4-4 

K  =  12’16 

11 

12 

4.5 

D 

4-5 

K  =  2-5 

7 

8 

4.5 

4-6 

K  =  2"  9 

i! 

12 

4.5 

5.0 

4-9 

4-10 

With  reference  to  Figures  4-2  through  4-4,  the  frequency,  f0,  during 
steady  state  averaged  from  6  to  11  Hz.  The  lowest  average  frequency  was 
obtained  with  the  condition:  qr  =  q to  =  0.  However,  for  the  no  leakage  con¬ 
dition  the  peak  to  peak  amplitude  oty  approached  0. 07  degree  as  compared 
to  approximately  0. 03  degree  with  leakage.  The  variation  of  and  f0  pro¬ 
duced  with  the  assumed  worst  case  leakage,  did  not  noticeably  differ  from 
the  variation  obtained  with  the  hardware  parameters. 
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Figure  4-2.  Recorded  Transients  for  an  Initial  Yaw  Condition, 
yo  -  -1/2°  (Ref:  Table  4-1) 
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Figure  4-3.  Recorded  Transients  for  an  Initial  Yaw  Condition 
(i//Q  =  -1/2°)  and  no  Leakage  (qL  ^  qLC  -  0) 
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Figure  4-4.  Recorded  Transients  for  an  Initial  Yaw  Condition 
('I'o  =  -1/2°),  Assuming  the  "Worst”  Case  Leakage 
(qL  =  0.150  in’/s,  qT^,  =  -0. 0165  inVs) 


2.  EFFECT  OF  COMPENSATION  GAIN,  K 

System  performance  and  response  characteristics  of  the  missile  can 
be  readily  modified  by  changing  the  gain,  K,  of  the  autopilot  compensation. 
Transients  for  K  =  2-16  are  shown  in  Figure  4-5  and  for  K  =  2-5  in  Figure 
4-6.  As  K  approaches  zero  (K  =  2-16)  a  steady  state  error  of  ^  results  due 
to  the  presence  of  leakage.  For  large  K  (K  =  2-5)  the  missile  transient  is 
highly  oscillatory. 

3.  COMPARISON  WITH  ANALOG  FLIGHT  CONTROL 

Table  4-3  shows  the  performance  of  the  analog  flight  control  system 
simulated  as  part  of  the  Phase  1  study.  Transients  for  a  yaw  step  command 
and  an  initial  yaw  condition  are  shown  respectively  in  Figures  4-7  and  4-8. 
Similar  performance  and  transient  characteristics  are  achieved  with  the  digi¬ 
tal  flight  control  with  the  compensation  gain  K  equal  to  2-9  (Figures  4-9  and 
4-10).  The  maximum  stable  step  for  the  analog  system  is  0. 075  degree 
larger,  a  difference  of  8  percent.  This  is  attributed  to  the  change  in  imple¬ 
mentation,  (see  Table  2-1).  With  the  digital  flight  control,  oscillations  of 
V'  during  steady  state  approached  a  peak- to- peak  value  of  0. 04  degree,  a 
consequence  of  a  nonlinearities  and  quantization.  Because  of  the  dynamics 
of  the  vehicle,  the  effect  of  the  oscillations  of  ^  on  the  flight  path  angle  is 
attenuated. 

As  discussed  in  Reference  1,  the  no- input  condition  used  to  investigate 
gain  of  the  digital  flight  control  is  not  analogous  to  the  no- input  condition  used 
with  the  analog  flight  control.  A  0. 109  degree  oscillation  of  the  nozzle  exists 
during  the  steady- state  condition  of  the  digital  system.  Hence,  the  maximum 
gain  of  the  analog  system  is  higher.  The  gain  of  the  analog  system  approaches 
that  of  the  digital  system  when  measured  with  an  input  of  Ny  =  1 .  Although 
maximum  gain  of  the  digital  system  is  lower,  the  gain  changed  required  to 
produce  Instability  provides  a  credible  margin. 

4.  GAIN- PHASE 

The  frequency  response  and  phase  shift  characteristics  of  the  closed 
loop  system  is  depicted  in  Figure  4-11  when  using  the  parameter  values  of 
Table  4-1,  To  minimize  the  effect  of  velocity  saturation  the  input  was  re¬ 
stricted  to  small  signals. 
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Figure  4-5.  Recorded  Transients  for  a  Yaw  Step  Command 
(Ny  =  10)  Showing  Steady  State  Error  of  ^  when  K  —  0 

(K  =  2’16) 
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Figure  4-6.  Recorded  Transients  for  a  Yaw  Condition 
(fQ  =  -1/2°)  with  K  =  2-5 
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Table  4-3 


RESULTS  OF  ANALOG  FLIGHT  CONTROL 
SIMULATION  (PHASE  I) 


Run 

Input 

Results 

Time 

Histories 

C-13A 

C-13B 

None 

None 

Stable  with 

Kg  = 18. 2 

Unstable  with 
Kg- 18.  9 

C-14 

Step  of  Ny  =  12 

Stable 

C-15 

Step  of  Ny  =  13 

Unstable 

C-16 

Initial  =  0.  5° 

Stable 

Figure  4-8 

C-17 

Step  of  Ny  -  10 

Stable 

Figure  4-7 

C-18 

Step  of  Ny  =  1 

Stable 

C-19A 

Step  of  N  »  1 

Stable  with 

K  *7.0 

g 

C-19B 

Step  of  Ny  =  1 

Unstable  with 
Kg.  7. 25 

C-20A 

Step  of  Ny  =  2 

Stable  with 

K  =  4.  5 

K 

C-20B 

Step  of  N  =  2 

y 

Unstable  with 
Kg.  4. 75 

4-11 
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Figure  4-7.  Recorded  Transients  Using  Analog  Flight  Control 
for  a  Yaw  Step  Commanc  (Ny  =  10) 
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Figure  4-8.  Recorded  Transients  Using  Analog  Flight  Control 
System  for  an  Initial  Yaw  Conation  (\pQ  =  +1/2°) 
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Figure  4-9.  Recorded  Transients  for  a  Yaw  Step  Command 
(Ny  =  10)  with  K  =  2" 9  (Ref  -  Table  4-1) 
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>-10,  Recorded  Transients  for  an  Initial  Yaw  Condition 
W0  = -1/2°)  with  K  =  2-9  (Ref  -  Table  4-1) 
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Section  V 


LINEARIZED  MODEL 


As  a  first  approximation,  a  linearized  model  provides  a  convenient 
means  of  determining  the  suitability  of  the  DHA  to  a  given  application.  To 
ascertain  the  validity  of  using  a  linear  model,  the  mathematical  model  of 
Figure  2-3  is  reduced  to  a  linear  form.  Performance  of  the  linear  model  is 
compared  to  the  performance  of  the  simulated  system  when  closing  the  loop 
with  the  laboratory  model.  Validity  is  established  by  the  agreement  of  per¬ 
formance  characteristics  for  frequency  response,  phase  shift,  and  step  re¬ 
sponse. 

A.  DEVELOPMENT 

Figure  5- 1  shows  the  block  diagram  of  the  digital  flight  control  system 
when  load  dynamics  are  neglected.  Restricting  the  commands  so  that  the  in¬ 
put  to  the  non-linear  block  never  exceeds  one  quantum  and  neglecting  quan¬ 
tization  Figure  5-2  is  equivalent  to  Figure  5-1.  By  use  of  the  identity  shown 
in  Figure  5-3  the  model  can  be  represented  by  the  diagram  of  Figure  5-4 
which  also  includes  the  approximation: 


Ll±  -v 

TjS 


2 


e 


-TjS 
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Zero- order  holds  can  be  formed  in  the  feed  forward  path  and  minor 
loop  by  combining  the  integration  of  the  DHA  with  the  first  difference  compu¬ 
tation.  The  sampling  frequency  is  approximately  two  orders  of  magnitude 
higher  than  the  frequency  of  the  missile  poles.  Hence,  a  sample  and  hold  car 

-TS 

be  approximated  by  the  time  delay,  e  Using  this  approximation  and  re¬ 
arranging  the  terms  yields  the  linearized  model  shown  in  Figure  5-5. 

The  linearized  model  is  based  on  small  signal  inputs,  i.  e. ,  velocity 
saturation  does  not  exist.  The  magnitude  of  missile  poles  are  small  relative 
to  the  poles  of  the  actuator.  As  a  consequence  loop  dynamics  even  in  the 
presence  of  slew  saturation  are  dominated  by  the  characteristics  of  the  mis¬ 
sile.  Therefore,  the  linearized  model  can  be  used  with  inputs  which  produce 
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Figure  5-1.  Block  Diagram  of  Digital  Flight  Control  System 
Neglecting  Load  Dynamics 


Figure  5-2.  Block  Diagram  Equivalent  of  Figure  5-1  When 
Neglecting  Quantization  and  Restricting  Input 
to  Small  Signals 


Figure  5-5.  Linearized  Model  of  Digital  Flight  Control  System 


5-3 


velocity  saturation.  The  magnitude  of  the  input,  however,  is  limited  by  the 
stability  condition: 


[Ts  ■toty 

i  &  > 

|  max 

CN/3  *S  ‘p/3 

(t) 

and  |t2  sin  £  lfi  j 

8(t)  > 

CHP  'qS  V 

(t)  for  sgn 

B.  RESULTS 

8(t)  =  -  sgn  (3  (t) 


Results  presented  in  Figure  5-6  and  5-7  were  obtained  using  the  line¬ 
arized  model  of  Figure  5-  5.  Although  not  included  in  this  report,  satisfactory 
results  were  also  achieved  with  a  model  which  included  the  approximation: 


which  was  used  to  obtain  the  response  curve  of  Figure  5-7.  Comparison  of 
Figure  5-6  and  5-7  with  Figures  4-11  and  4-1  shows  the  linearized  model  to 
be  a  good  approximation. 


Figure  5-6.  Closed  Loop  Amplitude  and  Phase  Characteristics 
of  Linearized  Model  (Ref  Fig.  5-  5) 
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Figure  5-7.  Computed  Response  to  a  Yaw  Command  (N  =  10) 
Using  Linearized  Model  (Figure  5-5)  ^ 
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Appendix  A 

DIGITAL  HYDRAULIC  ACTUATOR 
OPERATIONAL  CONCEPT 


The  digital  hydraulic  actuator  (DHA)  design  concept  used  in  this  study 
is  shown  in  Figure  A-l.  Calibrated  fluid  volumes,  transferred  by  displacing 
the  calibrating  piston,  are  summed  to  achieve  a  load  position  corresponding 
to  a  digital  command.  Auxiliary  pistons  permit  the  calibrated  volume  to  be 
displaced  at  a  rate  independent  of  the  mass  load. 

For  simplicity,  the  operation  of  the  DHA  (Figure  A-2)  will  be  described 
neglecting  the  auxiliary  pistons.  The  high  pressure  is  directed  to  force  the 
volume  piston  to  its  extreme  upper  position.  When  the  actuating  valve  is 
moved  down,  the  fluid  pressure  will  move  the  volume  piston  to  its  extreme 
bottom  position  and  displace  fluid  into  the  actuating  cylinder.  This  in  turn 
will  cause  the  output  ram  to  move  a  calibrated  distance  to  the  right.  Because 
of  the  larger  area  on  the  left  side  of  the  ram  piston,  the  high-pressure  fluid 
on  the  right  side  of  the  cylinder  will  be  forced  into  the  high-pressure  line.  A 
restriction  placed  in  this  line  will  damp  the  output  ram.  Further  movement 
of  the  actuating  valve  will  cause  the  output  ram  to  continue  discrete  motions 
to  the  right. 

To  see  how  motion  to  the  left  is  achieved,  assume  that  both  valves  are 
again  in  the  upper  position,  but  this  time  move  the  direction  valve  down.  The 
return  pressure  will  now  propagate  to  the  bottom  side  of  the  volume  piston, 
causing  the  volume  piston  to  move  down.  This  will  withdraw  one  calibrated 
volume  from  the  actuating  cylinder,  and  the  output  ram  will  be  forced  to  the 
left  by  the  high  pressure  on  the  right  side  of  the  actuator  cylinder.  Motion 
to  the  left  is  continued  by  moving  the  actuating  valve.  Hence,  the  output  ram 
is  moved  to  the  right  by  moving  the  direction  valve  up,  or,  if  it  is  already  in 
its  upper  position,  by  changing  the  position  of  the  actuating  valve.  Output 
ram  motion  to  the  left  is  acl  ieved  by  moving  the  direction  valve  down,  or,  if 
it  is  already  down,  by  changing  the  position  of  the  actuating  valve. 

Accuracy  requires  that  the  valve  iteration  interval,  Tc,  oe  greater  than 
the  time  required  to  transfer  the  calibrated  volume,  Vc.  The  maximum  valve 
iteration  rate,  fc,  for  a  given  pressure  is  thus  dictated  by  the  magnitude  of  the 
mass  load.  The  addition  of  auxiliary  pistons  (Figure  A-l)  permits  the  cali¬ 
brated  volume  to  be  displaced  at  a  rate  independent  of  the  mass  load. 
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Figure  A-l.  Schematic  of  the  Basic  Digital  Hydraulic  Actuator 


Figure  A-2  Digital  Hydraulic  Actuator 


Hence,  the  iteration  rate  is  limited  only  by:  (1)  the  maximum  switching  speed 
of  the  valve;  (2)  the  time  required  to  transfer  Vc;  and  (3)  the  time  required  to 
decay  the  transient  pressures  displacing  the  volume  piston.  Since  the  mass 
of  each  auxiliary  piston  is  much  less  than  that  of  the  mass  load,  fluid  transfer 
time  is  greatly  reduced.  Consequently,  the  mass  rate,  mr,  acting  to  displace 
the  mass,  M,  approaches  a  step  function. 


If  the  actuating  valve  is  positioned  so  that  high  pressure,  Ph,  forces 
the  calibrating  piston  from  the  position  shown  in  Figure  A-l  to  the  bottom 
position,  a  calibrated  volume  of  fluid,  Vc,  will  be  displaced.  The  volume 
displaced  can  be  considered  as  being  divided  into  three  unequal  parts: 


•  A  volume  necessary  to  increase  the  transient  total  internal  pressure, 
PL,  from 


A  volume 
A  volume 


to  P.  a2  . 

*2 

which  displaces  the 
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ram  piston  the  right, 
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Because  the  mass,  M,  driven  by  the  ram  piston  is  large  relative  to  the 
mass  of  the  auxiliary  piston,  the  response  of  the  auxiliary  piston  is  high. 
Hence,  the  larger  portion  of  the  displaced  volume  is  used  to  displace  the  aux¬ 
iliary  piston.  The  ram  piston,  initially  accelerated  by  a  portion  of  the  volume 
displaced  by  the  calibrating  piston,  is  further  displaced  by  the  transfer  of  the 
volume  stored  by  the  auxiliary  piston  due  to  P.  acting  to  restore  auxiliary 
piston  2  to  its  rest  position.  Motion  of  the  ram  piston  is  completed  when  the 
volume  required  to  increase  P.  is  transferred  to  ram  piston  motion,  thus 
returning  Pl  to  its  original  vame,  P^  Ar/Al.  Hence,  the  ram  piston  dis¬ 
placement  Xj  is: 


Additional  calibrated  volumes  displaced  by  the  calibrating  piston  are  stored 
by  displacing  auxiliary  piston  2. 

Over -travel  of  the  ram  piston  as  it  moved  to  the  right  reduced  P.  .  As 
a  consequence,  auxiliary  piston  3  is  displaced,  thus  limiting  P.  to  minimize 
the  effect  of  cavitation.  Displacement  of  auxiliary  piston  3  creates  a  pressure 
differential  that  acts  to  return  the  auxiliary  piston  to  its  rest  position.  Hence, 
over -travel  of  the  ram  piston  is  limited.  Auxiliary  piston  3  is  restored  to  its 
rest  position  and  returns  the  ram  piston  to  a  position  such  that 


Positioning  the  sign  valv.  so  that  the  return  is  connected  to  the  cali¬ 
brated  volume  permits  the  ram  piston  to  be  displaced  to  the  left.  Displacing 
the  actuating  valve  allows  Pl  to  displace  the  calibrating  piston.  Due  to  the 
magnitude  of  the  mass,  M,  the  Initial  volume  of  fluid  required  to  displace  the 
calibrating  piston  is  furnished  by  the  motion  of  auxiliary  piston  3.  Operation 
is  therefore  identical  in  principle  to  the  motion  of  the  ram  piston  to  the  right. 
However,  when  ram  piston  motion  is  to  the  left,  auxiliary  piston  3  provides 
the  storage  capability  and  auxiliary  piston  2  returns  the  ram  piston  to  its 
commanded  position  when  over -travel  occurs. 

Hydraulically  actuated  pilot  valves  operate  both  the  sign  and  actuating 
valves.  In  addition,  by  using  hydraulic  power  to  move  the  pilot  valve,  the 
magnet  is  not  required  to  perform  any  work,  and  thif  permits  the  use  of  small, 
low-power  magnets  controlled  by  electrical  commands.  Energized,  the  mag¬ 
nets  act  as  latches,  holding  the  spool  of  the  pilot  valve  at  a  fixed  position. 

When  de-energized,  a  magnetic  latch  permits  the  hydraulic  force  to  move  the 


A-  Ti 


pilot  spool,  which  in  tu-n  allows  a  hydraulic  force  to  move  the  main  spool. 

A  schematic  diagram  ol  the  pilot  and  actuator  valves  is  shown  in  Figure  A-3. 
The  pilot  valve  is  shown  held  in  its  left  hand  position  by  magnetic  latch  A. 
Hydraulic  pressure  acts  through  the  pilot  valve  to  the  end  of  the  actuator  valve, 
forcing  it  to  the  extreme  right.  It  also  acts  through  the  actuator  valve  to  the 
left  end  of  the  pilot  spool.  De-energizing  magnetic  latch  A  allows  the  hydrau¬ 
lic  pressure  to  move  the  pilot  to  the  extreme  right-hand  position,  where  it  is 
held  by  latch  B.  Moving  the  pilot  to  the  right-hand  position  switches  the  pres¬ 
sure  to  the  right-hand  end  of  the  actuator  valve  and  connects  the  left  end  to 
the  sump.  The  actuator  valve  is  therefore  moved  to  the  extreme  left.  In  an 
extreme  left  position,  the  pressure  switches  to  force  the  pilot  to  the  left. 

The  left  end  of  the  pilot  is  now  connected  to  the  sump,  and  movement  of  the 
pilot  is  prevented  by  magnetic  latch  B.  De-energizing  latch  B  allows  both 
valves  to  move  back  to  the  initial  position. 
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Figure  A-3.  Pilot  Drive  Schematic 
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Description  of  Symbols 


Symbol 

Description 

Type  -  IBM  No. 

H 

□ 

h 

Inverter 

370348 

1 

□ 

h 

2-Input 

NAND 

Inverter 

370216 

ILL 

— 

h 

3-Input 

NAND 

Inverter 

370378 

H 

□ 

P 

Trigger 

37022 

H 

a 

h 

Singleshot 

370262 

H 

□ 

pr 

Driver 

371633 

hr 

Load  Card 

37032 

B-2 


Time  Relation 
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Figure  B-l.  Clock  Pulses  -  Function  and  Time  Relation  (Revised) 
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Figure  B-2.  Central  Logic 
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Figure  B- 14.  Large  Actuating  Driver 
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Symbol 

Dimension 

Definition 

Ac 

in2 

Area  of  Calibrating  Piston 

al 

in2 

Area  of  Ram  Piston  Exposed  to  P^ 

in2 

Area  of  Ram  Piston  -  Rod  End 

A2 

in2 

Area  of  Auxiliary  Piston  2  Exposed  to  PL 

a2 

in2 

Area  of  Auxiliary  Piston  2  Exposed  to  Pg 

A3 

in2 

Area  of  Auxiliary  Piston  3  Exposed  to  PL 

a3 

.  2 
in 

Area  of  Auxiliary  Piston  3  Exposed  to  Pg 

C 

lb  -  s/in  1 

Viscous  Damping  Coefficient 

CN/3 

— 

Gradient  of  Normal  Force  Coefficient  with 
Respect  to  Angle  of  Sideslip 

eb 

deg 

Error  Required  to  Produce  Input  -  Large 

es 

deg 

Error  Required  to  Produce  Small  Volume 
Command 

FB 

— 

Feedback 

Af 

— 

External  Force 

K 

— 

Forward  Path  Gain  of  Integration  Compensa¬ 
tion  -  Equal  to  K.T 
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uuu 


1 

fl 

3 
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Symbol 

Dimension 

Definition 

Ka 

— 

Pitch  and  Yaw  Attitude  Autopilot  Gain 

*b 

- 

Weighting  Factor  for  Yaw  Body  Angles 

Kg 

— 

Variable  Gain  Used  to  Determine  Gain  Margin 

Ki 

-1 

s 

Forward  Path  Gain  of  Autopilot  Integration 
Compensation 

kl 

deg  s” 1 

Leakage  Rate  Gain 

Kr 

s 

Yaw  Rate  Autopilot  Gain 

K2 

deg/pulse 

Proportionality  Factor  for  Guidance  Yaw 
Command  and  Pulses 

ft 

Aerodynamic  Control  Arm 

ft 

Axial  Control  Arm 

M 

lb  -  s2/in" 1 

Lumped  Mass  of  Output  Ram  Piston  Including 
Load  Inertia 

My 

radians/ 

pound 

Missile  Transfer  Function  i///Ty 

Ny 

— 

Algebraic  Sum  of  Yaw  Steering  Pulses 

NyFB 

kb 

TP5-* 

N 

y« 

N  (t)  +  N 
y  yFB 

pd 

lb/in2 

Pressure  Change  of  P. 

La 

ph 

lb/in2 

High  Pressure 

PL 

lb/in2 

Transient  Total  Internal  Pressure 
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Symbol 

Dimension 

Definition 

Pm 

lb/in2 

Maximum  Deviation  of  PL 

Pr 

lb/in2 

Return  Pressure 

Ps 

lb/in2 

Supply  Pressure 

q 

in3  -  s” 1 

Flow 

Q 

^corr 

in3/s 

Correction  Required  to  Make  Analog  Integrated 
Value  Vg  Equal  to  Integrated  Value,  Vdig 

in3/s 

Flow  for  Specified  Parts  of  the  Period  T 

Q 

P 

in3/s 

Summation  Q_.  and  Q  at  Sample  Time 

PL  corr 

q 

lbf/in2 

Dynamic  Pressure 

in3/s 

2 

[r  Vprl 

qLC* 

“LC  * 

|h^7  — 5— J 

qLC 

inV1 

Error  Rate  Leakage 

qL 

in3/ s” * 

Leakage  Rate  Gain 

rl 

lb  s/in"3 

Flow  Resistance 

R 

in 

Moment  Arm 

S 

s"1 

Laplace  Transform  Variable 

7 

ft2 

Missile  Reference  Area 

T 

s 

Sampling  Period 

td 

s 

Delay 
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Symbol 


TS 

AV 


AV' 

Vdig 

AVi 

AV2 


Dimension  Definition 

s  Valve  Operation  Period 

s  Fluid  Flow  Period  (Large  Volume) 

s  Fluid  Flow  Period  (Small  Volume) 

lb  Thrust  Along  Yaw  Axis 

lb  Nominal  Total  Thrust 


in 


.  3 

in 


in 


in 


in 


in 


in 


in 


.  3 
in 


Volume  Transferred  Over  Interval  t  to 
(tv  +  T) 

Volume  Transferred  Over  the  Interval  t . 
t0  ^A/D+  T>  A/D 

Digital  Summation  of  AV' 

Volume  Transferred  Over  the  Interval  t 
to  (tA/D  ♦  T) 

Volume  Transferred  Over  the  Interval 
«A/D  *  T)  to  (ty  +  T) 

Analog  Integrated  Value  of  the  Flow 

Leakage  Volume  Change  Due  to  the  Pressure  P^ 

Total  Volume  Producing  Ram  Motion 

Volume  Generating  Pressure  Pd 
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Symbol 

Dimension 

Definition 

VCS 

in3 

Displacement  Volume  of  Calibrating  Small 
Piston 

V 

in3 

Oil  Volume  of  Digital  Hydraulic  Actuator 

VCB 

in3 

Displacement  Volume  of  Calibrating  Piston  - 
Large 

VT 

in3 

Total  Volume  Producing  Ram  Motion 

X1 

in 

Displacement  of  Output  Ram  Piston 

X1A 

in 

Actual  Displacement  of  Output  Ram  Piston 

8a 

deg 

Actual  Nozzle  Rotation 

h/3 

deg 

Feedback  Gain 

8c 

rad 

Desired  Nozzle  Position 

8n 

deg 

Feedback  Gain 

8 

deg 

Feedback  Gain 

P 

8L 

deg 

Feedback  Gain 

8lc 

deg  s  " 1 

Error  Rate  Leakage 

$ 

lbf/in2 

Bulk  Modulus 

m 

rad 

Angle  of  Sideslip 

\ 

deg 

One  Quantum 

s 

deg,  rad 

Rotatable  Nozzle  Cant  Angle 

+ 

deg,  rad 

Yaw  Angle 

f 

deg 

Desired  Change  of  Nozzle  Position 
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